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The quantum-classical crossover from the Fermi liquid towards the Wigner solid is numerically 
revisited, considering small square lattice models where electrons interact via a Coulomb U/r 
potential. The studies of models without disorder and spin and including disorder and spin 
show that the electron solid is formed in two stages, giving rise to an intriguing solid-liquid 
regime at intermediate couplings. 
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1 Lattice model 

We consider A^ electrons on L x L square lattice with periodic boundary conditions (BCs). The 
Hamiltonian reads 

W = ^(-t ^ 4,aCi,^ + Vini,a) + ^ Zl ':"_'':"' + 2f/ Z rii^^m^i, (1) 

^^ ■ i^i' 

a, 

Qj ' where q^o- (c|g.) destroys (creates) an electron of spin a at the site i and nj^o- = So-Ci.o-- The 

first terms describe the kinetic energy (oc —t) and the random substrate energy (potentials Vi 
uniformly distributed inside [—W/2,W/2\). The interaction consists of a U /\i — i'\ Coulomb 

"^ ' repulsion plus a 2U Hubbard repulsion. |i — «'| is the smallest distance between the sites i and 

^—N ■ z' on a square lattice with periodic BCs. In our model, the Coulomb energy to kinetic energy 

lO . ratio r^ = U/{2ty/Tm^) for a filling factor Ue = N/L?' . 

trr* \ S and Sz are the total spin and its component along an arbitrary direction z. Ti can be 

written in a block-diagonal form, with A^ -|- 1 blocks where Sz = —N/2, . . . , N/2 respectively. 
Assuming A^ = 4 and L = 6, the three blocks with ^^ > are diagonalized using Lanczos 
algorithm to obtain the minimum eigenenergy Eo(Sz) of each block. 

2 Spinless fermions without disorder 

We begin by studyingEJ the ground state (GS) of the block TC{Sz = 2) when W = 0. When [/ = 0, 
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the states are A'^^- plane wave Slater determinants (SDs) np=i^I('„)|0 >> where dL \ creates a 
particle in a state of momentum k{p) = 2Tr{px,Py)/L and |0 > is the vacuum state. The low 
^ ■ energy eigenstates are given by: (A) 4 degenerate ground states (GSs) \Ko{(3) > (/3 = 1, . . . ,4) 

H I of energy Eq(U = 0) = —13t and of momenta Kq ^ 0; (B) 25 first excitations of energy 

Ei{U = 0) = -I2t. (C) 64 second excitations 1^2 (a) > of energy E2{U = 0) = -lit; (D) 180 
third excitations |i^3(a) > of energy E^{U = 0) = —lOt and (E) 384 fourth excitations of energy 
E4{U = 0) = —9t. We define the 20 plane wave SDs useful to partly describe the intermediate 
GS. They are given by 4 plane wave SDs |i^i(/3) > of energy — 12t where the particles have 
energies —At, — 3t, —3t, —2t respectively, and by 16 plane waves SDs \K4^{6) > of energy —9t 
where the particles have energies —At, — 3t, —2t, Ot (first set of 8 SDs) and — 3t, — 3t, —2t, —t 
(second set of 8 other SDs) respectively. Those 20 SDs are directly coupled by the pairwise 
interaction and their momenta are zero {K = X]p=i k{p) = 0). 

When t = 0, the states are A'^^ Slater determinants cjcjc^cj |0 > built out from the site 
orbitals. The low energy states are the following site SDs: (A) 9 squares |5'o(/) > (/ = 1, . . . , 9) 



of side a = 3 and of energy Eo{t = 0) ss 1.80U; (B) 36 parallelograms |5'i(I) > of sides (3, vlO) 
and of energy ^ 1.85C/; (C) 36 other parallelograms \S2{I) > of sides {VW, VW) and of energy 
PS 1.97U and (D) 144 deformed squares |5'3(I) > obtained by moving a single site of a square 
|5/ > by one lattice spacing and of energy « 2U. 
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Figure 1: UPPER LEFT: Low energy part of the spectrum exhibiting a GS level crossing at rf . Inset: two 
first level spacings Ai/t (dashed) and A2/t (dotted) which become equal at rY and the perturbative result 
Ai/f = A2/t ^ 10392/rf valid when r^ -» oo (dot-dashed). UPPER RIGHT: GS projections Po{rs) onto a few 
plane wave SDs, given by the 4\Ko{l3) > (empty circle), the 4 \Ki{l3) > (empty square), the 64 \Ki{a) > (filled 
diamond), the 180 \K2ia) > (x), the 16 \K4,{d) > (asterisk) respectively. LOWER LEFT: GS projection Poo{rs) 
onto a few site SDs, given by the 9 squares 150(1) > (filled square), the 36 parallelograms |Si(/) > (asterisk), 
the 36 parallelograms |52(/) > (diamond), and the 144 deformed squares 15*3(7) > (left triangle) respectively. 
LOWER RIGHT: GS projection Poire) (asterisk) and P^{rs) (empty square) and total GS projection P (filled 
circle) onto the re-orthonormalized basis using the low energy eigenvectors of the two limiting bases. 



For the first low energy states, the crossover from the U = eigenbasis towards the t = 
eigenbasis is shown in Fig. ffl (upper left) when one increases the ratio r^. If we follow the 4 GSs 
Eo{rs = 0) {Kq 7^ 0), one can see a first level crossing at rf ps 9.3 with a non degenerate state 
(Kq = 0) which becomes the GS above rf , followed by two other crossings with two other sets of 
4 states with Kj 7^ 0. When r^ is large, 9 states coming from Ei{rs = 0) have a smaller energy 
than the 4 states coming from £'0(^5 = 0). The degeneracies ordered by increasing energy become 
(1,4,4,4, . . .) instead of (4,25,64, . . .) for r^ = 0. These 9 low energy states give the 9 square 
molecules |5'o(/) > when r^ — > 00. When r~^ is very small, the first 9 states correspond to a 
solid molecule free to move on a restricted 3x3 lattice, with an effective hopping term T oc tr~^. 
This gives 9 states of kinetic energy given by —2T{cos Ki{I) + cos Kt{I)) with Ki{I) = lirpifi 
and Kt{I) = 27rpt/3 {pi^t = li2,3). This structure with degeneracies 1,4,4 respectively and two 
equal energy spacings Ai and A2 appears (inset of Fig. || upper left) when r^ is larger than the 
crystallization threshold r^ ~ 28. The two characteristic thresholds r^ (level crossing) and rf 
(9 first states having the structure of the spectrum of a single solid molecule free to move on a 



3x3 square lattice) can also be detected by other methods given in Ref. [1]. 

To understand further the nature of the intermediate GS between rf and r^ , we have 
projected the GS wave functions l^o('^s) > over the low energy eigenvectors of the two eigenbases 
valid for U/t = and for t/U = respectively. Let us begin with the U = eigenbasis. Below 
rf , each of the 4 GSs |^o(rs) > with Kq ^ has still a large projection Po(^siO) = X]/3=i I < 
^0 (rs)|i^o(/5) > P over the 4 non interacting GSs. There is no projection over the 25 first 
excitations and smaller projections P^ij's^'l) and Po(^s)3) over the 64 second and 180 third 
excitations of the non interacting system. Above rf , the non degenerate GS with Kq = has a 
large projection 

4 

Po(r.,l) = El<^o(rs)|i^i(/3)>P (2) 

13=1 

which is equally distributed over the 4 excitations | Ki {(3) > of momentum Ki = and a second 

significant contribution 

16 

Po(r-.,4) =Y.\< ^o(r.)|i^4(<J) > |' (3) 

<5=1 

given by its projection onto the 16 mentionned plane wave SDs \K4{5) > of energy — 9t. Above 
rf , its projections onto the 4 \Kq{P) >, the 21 other first excitations and the second and 
third excitations of the non interacting system are zero or extremely negligible. The total GS 
projection Po(rs) = Po(^s,l) + Poi^sA) onto the 4 |i^i(/3) > and 16 \K4{6) > is given in Fig. |l| 
(lower right) when r^ > rf . This shows us that a large part of the system remains an excited 
liquid above rf . 

We now study the GS projections Poo onto the t = eigenbasis. The GS projection 

9 

Poo{rs,0) = J2\ <nirs)\SQ{I) > \^ (4) 

7=1 

onto the 9 square site SDs \So{I) > is given in Fig. || (lower left), together with the GS projection 
PooifsjJ) onto the site SDs corresponding to the J degenerate low energy excitations of the 
t = system. The total GS projection P4)(''s) = J2p=oPoo{rs,p) onto the 9 squares |5'o(/) >, 
the 36 parallelograms \Si{I) >, the 36 other parallelograms 1 5*2(1) > and the 144 deformed 
squares |5'3(/) > is given in Fig. 11] (lower right) when r^ > rf. This shows us that the ground 
state begins to be a floppy solid also above rf . 

The site SDs and plane wave SDs are not orthonormal. After re-orthonormalization, the 
total projection P of |^o(''s) > over the subspace spanned by the 4 \Ki{(3) > and 16 |i^4((5) > 
and 225 site SDs of lower electrostatic energies are given in Fig. |l| (lower right), showing that 
|^o('"s) > is almost entirely located inside this very small part of a huge Hilbert space for 
intermediate r^, spanned by low energy SDs of different nature, and adapted to describe a solid 
entangled with an excited liquid. 

3 Magnetization and polarization energies in presence of a random substrate 

We now consider weakly disordered samples when the spin degrees of freedom are included a. 
Their role and the consequences of an applied parallel magnetic field which aligns only the spins 
without inducing orbital effects, have been the subject of Ref. [2]. The study of a statistical 
ensemble of samples with W = 5, N = 4 and L = 6 provides complementary signatures of a 
particular intermediate behavior. Let us note M the fraction of clusters with S* = 1 at P = 0, 
Q2 = Eo{Sz = 2) — Eq{Sz = 0) and Qi = Eq{Sz = 1) — Eq{Sz = 0) the Zeeman energies 
necessary to yield S = 2 and 5 = 1 respectively for a cluster with 5 = 0. 
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Figure 2: As a function of Vs, fraction M of clusters with S = 1 at S = (filled diamond, right scale), partial 

< log Qi > (filled circle, left scale) and total < log Q2 > (empty triangle, left scale) energies required to polarize 

S = clusters to S = 1 and 5 = 2 respectively. The straight line corresponds to 0.25 — 21ogr3. 



In Fig. ^, M is given as a function of r^. One can see a first threshold at r^ ~ 0.35 where 
the interaction can drive 5 = 1 in certain samples. Above a second threshold rf , M regularly 
decreases to reach a zero value at a third threshold r^^ ~ 9 where an antiferromagnetic square 
molecule is formed. The ensemble averages < logQi > and < log (^2 > (without taking into 
account the 5 = 1 spontaneously magnetized clusters) define the typical fields B necessary to 
yield S' = lor5' = 2ina5 = cluster. In Fig. 0, one can see an intermediate regime again for 



r[^ < Ts 
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where < logQi > becomes roughly independent of Ts, while Q2 oc r^ 



4 Conclusion 

One concludes that mesoscopic Wigner crystallization proceeds in two stages. In a clean system, 
a minimal description of the intermediate GS requires to combine the low energy states of the 
two limiting eigenbases. In this sense, the intermediate GS is neither solid, nor liquid, but 
rather the quantum superposition of those two states of matter. This is strongly reminiscent of 
the conjecture proposed by Aridreev and Lifshitzd for the quantum melting of a solid. A path 
integral Monte Carlo approachs has recently shown that a few electrons confined in a harmonic 
trap crystallize also in two stages, firstly via a radial ordering of electrons on shells and secondly 
via the freezing of the intershell rotation. As one can see, a two stage crystallization is not only 
characteristic of a mesoscopic harmonic trap, but also occurs in a mesoscopic 2d torus. To add 
a random substrate defavors the liquid state. The magnetization gives also the signature of an 
intermediate regime where Stoner ferromagnetism is defavored by Wigner antiferromagnetism. 
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